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Abstract 

We give a detailed general description of a recent geometrical dis- 
cretisation scheme and illustrate, by explicit numerical calculation, the 
scheme's ability to capture topological features. The scheme is applied 
to the Abelian Chern-Simons theory and leads, after a necessary field 
doubling, to an expression for the discrete partition function in terms 
of untwisted Reidemeister torsion and of various triangulation dependent 
factors. The discrete partition function is evaluated computationally for 
various triangulations of S 3 and of lens spaces. The results confirm that 
the discretisation scheme is triangulation independent and coincides with 
the continuum partition function. 

1 INTRODUCTION 

A very useful way to regularize a quantum field theory is provided by the lat- 
tice formulation introduced by Wilson M. (See for example, || for a detailed 
treatment). However, this formulation has difficulty in capturing topological 
features of a field theory, for example, the topological theta-term in QCD. It 
is therefore of interest to investigate alternative discretisation schemes. In this 
paper we describe an alternative scheme which is applicable to antisymmetric 
tensor field theories including Abelian gauge theories and fermion field theory in 
the Kahler-Dirac framework || , and which is well-suited for capturing the topo- 
logical features of such theories. It is based on developing analogies between the 
different types of fields and the way they appear in a quantum field theory with 
a corresponding list of discrete variables and operators. The method is valid for 
any arbitrary compact 3-manifold without boundary. We illustrate the scheme 
by applying it to the pure Abelian Chern-Simons gauge theory in 3 dimensions, 
and to a doubled version, the so-called Abelian BF gauge theory. In the latter 
case the topological features of the theory are completely reproduced by our 
discretisation scheme, even before taking the continuum limit. We illustrate 
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this by explicit numerical calculations of the discrete partition function when 
the spacetime is S 3 or is a lens space, L(p,l). 

The discretisation scheme involves several ingredients which are not used in 
the standard Wilson lattice formulation. These include a triangulation of the 
spacetime (i.e. a decomposition into hyper-tetrahedra rather then hyper-cubes) 
and a mathematical tool called the Whitney map 0. These have previously 
been used to discretise field theories in H [| where various convergence results 
(e.g. convergence of the discrete action to the continuum action) were estab- 
lished. Triangulations of spacetime are also used in discretisations of other 
quantum field theories, for example in quantum gravity (see for a review). 
The thrust of the present work is quite different. We set up the discretisation 
in such a way that the geometric structures of the continuum field theory are 
mirrored by analogous structures in the discrete formulation. As we will see, 
this requires a certain doubling of fields. With this doubling the topological 
features of Abelian Chern-Simons theory are completely captured by the for- 
mulation. An alternate scheme of discretisation proposed previously in Q which 
does not involve the doubling of fields fails to capture the topological features of 
the Abelian Chern-Simons theory as we will show numerically in Sectionp. A 
mathematical treatment of this discretisation scheme has been given in [g~ llj . 
A short version of some of the results has been published in |j| . 

Our aims in this paper are, firstly, to make the techniques and results of 
[|[ accessible to a wider audience, and secondly, to demonstrate a practical 
numerical implementation of the discretisation scheme. Numerical implemen- 
tations are developed for Abelian Chern-Simons theory defined on the three 
dimensional sphere, S 3 , and on the lens spaces L(p, 1) for p = 1, 2 and 3. Our 
discretisation scheme is the only one which numerically reproduces the exact 
topological results for Abelian Chern-Simons theory p"2[ . 

This paper is organised as follows. In Section || the discretisation scheme is 
described together with a summary of the topological results used to set it up. 
In Section |3| some features of the Abelian Chern-Simons theory on general 3- 
manifolds are reviewed. In Section ^ the discretisation scheme is applied to this 
theory. An expression for the partition function of the resulting discrete theory 
is derived in terms of the data specifying the triangulation of the spacetime. In 
Section [s] the numerical evaluations of the partition functions corresponding to 
these triangulations are presented. In Section [s] we summarise our conclusions. 



2 THE DISCRETISATION SCHEME 

Let us start by considering a general quantum field theory defined on an arbi- 
trary manifold M of dimension D. Suppose the theory has fields <j) p (x) where 
x £ M , and where <p p is a p-form (antisymmetric tensor field of degree p de- 
fined on M). The Lagrangian for the system involves the fields, the Laplacian 
operator, and possibly (as for the Chern-Simon theory) an antisymmetrized 
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first-order differential operator. 

In differential geometric terms, the theory is constructed using the following 
objects which are defined on the manifold M: p-forms </> p which are generalised 
antisymmetric tensor fields, the exterior derivative d : <j) p <fi p+1 , the Hodge 
star operator * : <f) p — > (f> D ~ p , which is required to define scalar products, and 
the wedge operator <jp A <fi q — <fi p+q . 

We want to construct discrete analogues of these objects. We begin by 
summarizing the basic properties of our operators of interest [jl3| . On a manifold 
M of dimension D, the operations (A, *, d) on p-forms, <fi p (p — 0, . . . , D), satisfy 
the following: 

1. 4> p A<j) q = {-l) pq 4> q A4> p . 

2. d((j) p A (j) q ) = d(j) p A <p q + {-l) p <j) p A d(j) q . 

3. *0P = 4> D - p . 

4. ** = (-l) Dp+1 . 

5. d 2 = 0, (d*) 2 = 0. 

6. d* = * d*, (d* is the adjoint of d). 
The following definitions will also be required 

• The Laplacian on p-forms A p = d p ^\d* + d* +1 d p . 

• The inner product (<p p , 4> p ) — J M (j) p A *<p p ' . 

A few examples might now be helpful. Firstly consider QED in 4 dimensions. 
The gauge field A — (f) 1 is a 1-form. The electromagnetic field is a 2-form given 
by F — dA. The action for the gauge field in QED is given by 

S(A) = (F, F) = {dA, dA) = I dA A *dA. (1) 

Jm 

Thus S involves the operators *, d and the wedge product. Similarly for Abelian 
Chern-Simons theory the gauge field A and electromagnetic field F = dA are 
1-forms and 2-forms respectively, as for QED. The action for the theory is given 

by 

S{A) = I A A dA = (A, *dA) (2) 
Jm 

where the spacetime M is 3-dimensional. Note that in both cases we can think 
of the action S(A) as a quadratic functional of the gauge field A. 

We would like to discretise the fields of <fi p , the inner product (•,•), and the 
operators (A, *, d) such that discrete analogues of their continuum interrelation- 
ships hold. To do this it is necessary to first introduce a few basic ideas of 
discretisation. We start by discretising the manifold M. This involves replacing 
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Figure 1: A 3-simplex. 



M by a collection of discrete objects, known as simplices, glued together. We 



Firstly, for p > 0, a p-simplex a^ p > = [vq, . . . ,v p ) is denned to be the convex 
hull in some Euclidean space M. D of a set of p+1 points «o, V\, . . . , v p £ M. m . Here 
the vertices v t are required to span a p-dimensional space. This requirement will 
hold so long as the equations, Y^=o ^i v i = an d J2i=o — admit only the 
trivial solution Xi = for i = 0, .. . ,p for real. 

A few examples might clarify the geometry. Consider = [vo\. This 
is a point or 0-simplex. Next = [vq,vi] is a line segment or 1-simplex. 
An orientation can be assigned by the ordering of the vertices, in which case 
— er^ 1 ) = [ui,Wo] for example. The faces of a 1-simplex are its vertices [vo] 
and [v\] which are 0-simplices. cr*- 2 -* — [1)0,1)1,1)2] is a triangle or 2-simplex. 
We note that an even permutation of the vertices has the same orientation 
as er^ 2 ) while an odd permutation reverses it and will be written as — rr' 2 '. 
The faces of a 2-simplex are its edges [i>ojKi]i [fijfa], and [f2, w o]- Finally, 
cr( 3 ) = [v ,vi,V2,v 3 ] is a tetrahedron or 3-simplex. Its faces are the four triangles 
[vo 1 v 1 ,v 2 ],[vo,V2,v 3 ],[vo,vi,v 3 ] and [vi,v 2 ,v 3 ] which bound it. 

An important feature of our discretisation scheme is that the original sim- 
plices are subdivided by using simplex barycentres. Geometrically the barycen- 
tre of a p-simplex, a^ p \ is the point which represents its "center of mass". We 
denote the barycentre of crW = [vo, ■ ■ ■ , v p ] as the point 



need a few definitions Q . 




(3) 



As an example, the barycentre of <r^ 



[vo, vi] is the midpoint of the line 
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segment which joins the vertices vo and v\ . 

We can now describe a particular way that a given manifold M can be 
discretised. Let S be a collection of simplices {cr^}, n = 0, 1, . . . , D, with the 
property that the faces of the simplices which belong to S also belong to it. 
The elements of S glued together in the following way is known as a simplicial 
complex 0,0]: 

1. erf™-' n = if cr-™\ <ffl have no common face. 

2. cr|™' > nerj fc ' ) ^ if <7, J -"' ) , er^ have precisely one face in common, along which 
they are glued together. 

In many cases of interest (including all 3-manifolds and all differentiable mani- 
folds 0|), M can be replaced by a complex K which it is topologically equiv- 
alent to. K is then said to be a triangulation of M(Note this triangulation is 
not unique). In this way of discretising M, the building blocks are zero, one, 
. . ., D-dimensional objects, all of which are simplices e.g. generalised oriented 
tetrahedra. 

We now observe that the same manifold can be discretised in many different 
ways. In the discretisation described, we used simplices. We could just as well 
have used generalised oriented cubes. 

There is another method of discretising a manifold which is the dual of the 
simplicial discretisation just described. It associates with a simplicial complex 
K a dual complex K. We proceed to describe this construction. We will see 
that the basic objects of the dual complex K are again zero, one, two, . . ., D 
dimensional objects, but this time they are no longer simplices. We illustrate 
the method by considering a manifold, M, which is a disc. This is a manifold 
with a boundary. We triangulate this by the simplicial complex, K, shown in 
Figure ||. 

Now consider the barycentres of the building blocks of the simplicial complex 
K. We have the following list shown in Table [l]. 

Pictorially the simplicial complex K with its barycentres is shown in Fig- 
ure ||. 

We can now construct the dual description of a triangulation K. Geomet- 
rically this utilises a discrete analogue of the Hodge * operator. Recall the * 
operator maps a p-form to a [D— p)-form , where D is the dimension of the man- 
ifold M on which the p-form is defined. In the dual geometrical decomposition 
of the manifold, we want to set up a correspondence between a p-dimensional 
object and a (D — p) dimensional object. 

This is done as follows. We first construct (D — p) dimensional objects 
whose vertices are barycentres of a sequence of successively higher dimensional 
simplices, where each simplex is a face of the following one. In other words 
(D — p) dimensional objects of the form {cr p , <7 p +i, . . . , <td}, where a n is a face 
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Figure 2: Triangulation K of a disc. 



Geometrical object in K dimension Corresponding barycentre 

Wi — n n r x°) 

(0) r 1 n ~(°) 

02 = M CT 2 = v 0, 

0f = N o i 4°) = , , 

01 1) = bo,^l] 1 = |(«0 + «l) = «3 

erf- 1 = W 2 ] 1 = + W 2 ) = V4 

^3 1} = [«2, Wo] 1 0i X) = §(«2 + V ) = V 5 

0i 2) = [vo,v 1 ,v 2 ] 2 ct£ 2) = |(u +v 1 +v 2 ) = v 6 

Table 1: Elements of the simplicial complex K and their barycentres. 



of a n +i- The orientation of these are set so as to be compatible with the 
manifold. Joining these objects together gives us the dual of a p . 
Thus for instance the map * acts on [v ] as follows. 

*k ■ M -> e i[v ,di,a] U [v , <f 3 , a]e 03 . 

Where the orientation of each of the small triangles has to be coherent with the 
orientation of the original triangle. This is shown in fig. 3 and leads to mapping 
[vo] to the shaded two dimensional region. The orientations of the simplices are 
specified by arrows in the figure. Coherence of orientation means, for example, 
that the arrow of an edge agrees with the arrow of the triangle to which it 
belongs. Next we consider [w ,Wi]. This is a 1-simplex and is to be mapped to 
a (2-l)=l dimensional object. The map is defined as 

*k : [vo,vi] — > [<fi,a]. 
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Figure 3: Barycentres of K. 



Again the orientation of [<fi , a] has to be coherent with the orientation of the 
triangles already introduced when the map for [v ] was considered. Similarly 

*K ■ [V1,V 2 ] -> [CT2,(T], 

K • [V2,V ] -> [cT3,a-], 

and finally 

*K ■ [fO,Ul)V2] -» [ff]- 

We then have the alternate discretisation K for M shown in fig. || 

Note when two edges are glued together they must have opposite orien- 
tations. We can now give the general rule for mapping an n-simplex a n = 
[vq, . . . ,v n ] to a (D — n) dimensional object ( (D-n) cell) as follows: We think 
of a n as an element of a simplicial complex K. We have 

*k ■ [vq, ...,v n ]—> U[a n , (Tn+l, C£»]j 

where a n+ i is the barycentre of an (n+l)-simplex which has face. (j n +2 

is the barycentre of an (n+2)-simplex which has <j n +i as a face and so on. These 
objects have to be coherently oriented with respect to [vq, .. . ,v n ]. The set of 
these cells constitutes the dual space K of K. 

By this procedure we claim, a discrete version of the Hodge star operation * 
has been constructed. Let us explain. The Hodge * operator involves forms. It 
maps p-forms in D dimensions to a (D-p)-form. The *k map involves not forms 
but geometrical objects. There is a simple correspondence relation between 
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Figure 4: Dual complex K. 



these two cases. Given a p-form, <p p , and a p-dimcnsional geometrical space, 
S p , the p-form can be integrated over S p to give a number. Thus S p and <f> p 
are objects that can be paired. We can write this as a pairing 




In order to proceed, we need to introduce some more structure. We start by 
associating with a simplicial complex K, containing {a l p } (i = 1, . . . ,K p ;p — 
0, . . . , D) a vector space consisting of finite linear combinations over the reals 
of the p-simpliccs it contains. This vector space is known as the space of p- 
chains,C p (.K'). For two elements <r p ,a p S C P (K), a scalar product ((J p ,a p ) = 6* 
can be introduced. An oriented p-simplex changes sign under a change of ori- 
entation i.e. if a p — [vo, ■ ■ ■ , v p ] and r is a permutation of the indices [0, . . . ,p], 
then [w T (o), • • ■ , v r(p)] = ( — l) r [ u 0; ■ ■ * , v p ], with r denoting the number of trans- 
positions needed to bring [v T (o), • • ■ , v r(p)] to the order [i> , . . . , v p ). 

Given the vector space C P (K), the boundary operator d K can be defined as 

d K : C P (K) -> C p -i(K). 

It is the linear operator which maps an oriented p simplex to the sum of its 
(p-1) faces with orientation induced by the orientation of <r p . If a v = [vq, . . . , v p ] , 
then 

p 

da? = J2(-mVo,-.-, vu... v P ], 
i=0 

where [v ot . . . , Vi, . . . , v p ] means that the vertex Vi has been omitted from a p to 
produce the face "opposite" to it. 
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Given that C P (K) is a vector space, it is possible to define a dual vector 
space C P (K), consisting of dual objects known as cochains; that is we can take 
an element of C P (K) and an element C P (K) to form a real number. Since the 
space C P (K) has a scalar product, namely if cr p ,a p £ C P (K) then (a p , a p ) = Sij. 
We can use the scalar product to identify C P {K) = C P (K), so that we can 
consider oriented p-simpliccs as elements of C P (K) as well as C p (K). We can 
write our boundary operation as 

([u , . . • ,v h . . . ,Vp] ) d K [v , . • • ,Vp\) = 

This suggests introducing the adjoint operation cIk defined as 

(d K [v , ...,Vi,.. .,v p ], [w , ■ • -,vp\) = ([v , ...,Vi,.. .,Vp],d K [v , . . .,v p ]). 

This is the coboundary operator which maps C P {K) — > Cp+\(K). 
Indeed we have 

d K [v , ...,v p )= ^2[v,v , . . .,v p ], 

V 

where the sum is over all vertices v such that [v, Vq, . . . , v p ] is a (p+1) simplex. 

The boundary operators &k and the coboundary operator have the prop- 
erty Bk9k = dftdx — 0. Furthermore, 

da ■ Cp — > Cp+ii 
9k ■ C p — > Cp-%. 

These operators are the discrete analogues of the operators d and (— 1)- d (p+ 1 )+ 1 * 
d* = d* which act on forms. 

These operators could be defined only when a scalar product ("metric") was 
introduced in the vector space C p 's. At this stage we have a discrete geometrical 
analogue of d, d* and *. We have also commented on the fact that the operation 
* maps simplices into dual cells i.e. not simplices. If the original simplicial 
system is described in terms of the union of the vector spaces of all p-chains 
then the space into which elements of the vector space are mapped by * is not 
contained within this space, unlike the situation for the Hodge star operation 
on forms. We will see that this difference leads inevitably to a doubling of the 
fields when discretisation, preserving topological structures, is attempted. 

We now need a way to relate a p-chain to a p-form. This together with a 
construction which linearly maps p-forms to p-simplices will allow us to translate 
expressions in continuum QFT to a corresponding discrete geometrical objects. 
We start with the construction of the linear maps from p-chains to p-forms due 
to Whitney @. 

In order to define this map, we need to introduce barycentric coordinates 
associated with a given p-simplex a p . Regarding a p as an element of some M. N , 
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we introduce a set of real numbers {jio, . . . , ji p ) with the property 

Mi > o, 

i 

A point x G a p can be written in terms of the vertices of a p and these real 
numbers as 

v 

X = 2J HiVi- 
i=Q 

Note if any set of fii = then the vector x lies on a face of <r p . One can think of x 
as the position of the center of mass of a collection of masses (/xo, . . . , fj, p ) located 
on the vertices (vq. . . . ,v p ) respectively. Setting fii = for instance means the 
center of mass will be in the face opposite the vertex Vi. The Whitney map can 
now be defined. We have 

W K : C P (K) -» $ P (A), 

where $p(A) is a p-form. If a p E C P (K) then 

p 

W[a p ] = p! ^(-1) A . . . dfii A . . . d^i pi 

i=0 

where dfii means this term is missing, and 
ordinate functions of a p . 

We next construct the linear map from 
as the de Rham map. We have 

A K : $ P (A) -» C P (K), 

defined by 

< A K (^ p ),a p >= [ $ p , 

JcrP 

for each oriented p-simplex G K. 

A discrete version of the wedge product can also be defined using the Whit- 
ney and de Rham maps such that A K : C P {K) x C q {K) -> C p+q {K) as follows: 

x A K y — A K {W K {x) A W K {y)). 

It has many of the properties of the continuous wedge product in that it is 
skewsymmetric and obeys the Leibniz rule but it is nonassociative. 

At this stage we have introduced all the building blocks necessary to discre- 
tise a system preserving geometrical structures. We summarize the properties 
of the maps introduced in the form of a theorem p] : 



(fj,o, . . . , fj, p ) are the barycentric co- 
p-forms to p-chains. This is known 
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Figure 5: Barycentric subdivision of K. 



1. A K W K = Identity. 

2. dW K = W K d K , where d : <p p -> (f> p+1 . 

3. f m W K {a) =<a,P>,a,PeK 

4. d K A K = A K d. 

This theorem shows how d K can be considered as the discrete analogue of d. 
We now show how * K can be considered as discrete analogue of *. For this we 
need barycentric subdivision. 

We recall that given a simplicial complex {of}, i=l,- • • , K p ;p = 0, . . . , D. A 
set of points(vertices) could be assigned to each simplex, namely erf. These are 
the barycentres. These vertices, regarded as vertices of a simplex, subdivide the 
original simpliccs to give a finer triangulation of the original manifold. This is 
a barycentric subdivision map BK. Clearly the procedure can be repeated to 
give finer and finer subdivisions in which the simplices become "smaller" . The 
procedure is illustrated for a triangle in Fig. 5. 

Note that all the barycentres are present as vertices of the barycentric subdi- 
vision and that * K acting on simplices belonging to the simplicial complex, K, 
associated with [vo,vi,V2] leads to objects which are not, in general, simplices 
of [vo,v\,V2] but belong to a different space K. However both K and K are 
contained in the barycentric subdivision B\v§,v\,vi\. This is a crucial observa- 
tion. In order to construct the star map, two geometrically distinct spaces were 
introduced. The original simplicial decomposition K with its associated set of 
p-chains C P (K) and the dual cell decomposition K with its associated set of 
p-chains C P (K). These spaces are distinct. However both belong to the first 
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barycentric subdivision of K . This allows the use the * operation if we think 
of K and K as elements of BK. 

We proceed as follows. Let BK and K denote the barycentric subdivision 
and dual triangulation, and 

* K : C P {K) -> C n - p (K). 

However C P (K) and C P (K) are both contained in C P (BK) as we have seen. 
Let 

w bk . CP ( BK ) _> 4> P (M), 
denote the Whitney map. Then we have for x € C p {K),y e C ,n-p_1 (A')(l||: 



2. <9 A ' = (-1)^+! * k d k * K on C p (iT). 



These are the discrete analogues of the interrelationships between d,d* , * and 
< •, • > in the continuum. 

Note K ^ K and that properties of &k, Ak analogous to those for differential 
forms only hold if K , K are both regarded as elements of BK. This feature of 
the discretisation method is, as we shall see, crucial if we want to preserve topo- 
logical properties of the original system. If a discretisation method is introduced 
without the * operation in it then as we shall see in section ^| the topological 
properties of the partition function for the Abelian Chern-Simons gauge theory 
do not hold. 

We proceed to apply these ideas to the Abelian Chern-Simons gauge the- 
ory on a compact three manifold M. First we summarize properties of the 
continuum field theory. 

3 Schwarz's topological field theory and the Ray- 
Singer torsion 

We begin our treatment of continuum field theory by describing Schwarz's 
method for evaluating the partition function of the Chern-Simons gauge the- 
ory on a three-manifold M |l6| . We assume that the first real homology group 
( to be defined shortly) of the manifold vanishes; this is done for the sake of 
simplicity. Schwarz's method is applicable for arbitrary compact 3-manifolds 
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without boundary. Such manifolds are called "homology 3-spheres" . The main 
example we have in mind are the 3-sphere S* 3 and the lens spaces L(p, 1), p=l, 
2, . . .( for a definition and the basic properties of lens spaces see | 
fields of the theory are the 1-forms on M i.e. u> £ J7 1 (M). (In terms of a local 
coordinate system on M we have lo(x) — u^(x)dx fJ - .) The action is 

S(u) = [ ui A d±ijj — f dx 1 dx 2 dx 3 e fnyp uj ll d iy uj p . (4) 
Jm Jm 

Here and in the following Q, q {M) denotes the space of q-forms on M (i.e. the 
antisymmetric tensor fields of degree q) and d q : Q, q {M) — > f2 g+1 (M) i.e. the 
exterior derivative. It has the property d q d q -\ = so Im(<i g _i) C Ker(d g ), 
where Im(d q -\) is the image of the operator d q -\ while Ker{d q ) is the null 
space of the operator d q and the cohomology spaces H q (M) arc defined by 

H q (M) = Ker(d,)/Im(d g _i). 

The H q {M) are Abelian groups which contain topological information about 
the manifold. The vanishing of H 1 (M), for instance, holds if the manifold is 
simply connected that is any loop in M can be smoothly deformed to any other 
loop in M |l4| . Note that fl° (M) is the space of functions on M and since do is 
the derivative Ker(do) consists of the constant functions i.e. 

H°(M) = Ker(d o )/0 = Ker(d ) = R. 

Our requirement on M that H 1 (M) = implies that 

Im(d ) = Ker(dx). 

A choice of metric on M determines an inner product in the spaces Vl q {M) and 
allows the action (^) to be written as 

S(u>) = A < oj, (*di)uj >, (5) 

where * is the Hodge star operator. ( See jl3) for background on this and other 
differential-geometric constructions.) In order to evaluate the partition function 
of this action by Schwarz's method requires the introduction of the resolvent for 
S(lo). The partition function is defined as 

Z(X) = N J dioe lS{ul \ 

The main problem in evaluating Z(X) is to properly deal with the zeroes of S{lo). 
These zero modes contain topological information regarding the manifold as the 
space of zero modes is given by Ker(di) and hence should not be discarded. 
Schwarz introduced an algebraic method (the resolvent method) for dealing 
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with this problem. Although it is only valid for S(o;)'s which are quadratic in 
uj, it can be used to analyse 5(w)'s constructed on arbitrary compact manifolds 
without boundary. For systems of this type Schwarz's method is an algebraic 
analogue of the problem of gauge fixing. The advantage of the resolvent method 
is that it can be easily extended to deal with the process of discretisation as we 
will show. 

The resolvent is defined to be the following chain of maps 

-» R fi°(M) -^ d ° Im(do) = Ker(di) -> Ker(S) -» 0. (6) 

These chain of maps form an exact sequence, that is, the image of a map is 
the kernal of the map which follows. With the help of the resolvent, Schwarz 
was able to show that the partition function for the theory was given by jl6) 

Z{\) = e-T^(^)-§dct'((* ( i 1 ) 2 )-3det , Kdo) 1 det(c^o)-*, (7) 

where t is a non-topological geometry dependent function and £ as shown in 
[ pT~| is given by 

C = dim#°(M) - dim-ff^M). 

l is part of a phase factor and thus the absolute value of the partition function 
Z(X) is a topological quantity. This will be our main concern 

For completeness we give a quick proof of this result, ignoring phase factors 
and constants. 

Introducing a metric in the space of cj's allows us to write 
Z{\) = N [ duje lS{ul \ 

JKordie(Kcrdi)^ 

= Vol(Kerdi) • (detd^di)-^JV. 

We proceed to rewrite Vol(Kerdi) using the exact sequence associated with 
Kerdi and the manifold M. This procedure gives an expression for the partition 
function Z containing information about the spaces Kerdi and (Kerdi)- 1 -. Sim- 
ply dropping Vol(Kerdi) leads to a loss of information. We have Vol (KerS)= Vol (Kerdi) 
= Vol(Imdo) by assumption (if H\(M) is non-trivial, this equation has to be 
modified §). Also, 



do |(Kcr<io) i: ( Kerd o) ± -> (Imd ) 

Vol(Imd ) = | det'do | Vol(Kerdo)- 1 . 

Note 

n° = Kerdo © (Kerdo)- 1 
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and 



4>o ■ 7~to — ► Ho 
Vol(H ) = | det 0o | Vol (Ho). 

where Tio represents the space of harmonic 0-forms. Note harmonic p-forms are 
solutions of (d*d + dd*)(j> p = 0. In this space the Hodge star operator is present 
and hence a scalar product and volume can be defined. The map 4>o introduced 
relates the space of harmonic 0-forms to the space of de Rham cohomology Hq. 
By a theorem of Hodge this space of harmonic p- forms is isomorphic to the space 
of H p . The space of de Rham cohomology does not have a metric and hence 
we define the volume in this space with the help of the map <po- Therefore, 

Vol(Kerdo) 1 - = Vol(tto)[Vol(Kerdo)n\ 
= Vol(fto)(Vol(ff )r\ 
= (Vol(r!o))(det0o)~ 1 (Vol(Wo))" 1 . 

So that finally we get 

Vol(KerS') =| det'd || det0 o T 1 Vol(O )(VolH,) _1 . 

Choosing 

NVo\{Q. ){YoYH,)- 1 = 1, 

we get 

Z = Vol(KerS')(detdidi)"' =| det'd || Aetdxd* \ det^ I" 1 ■ 

A more careful calculation gives the determinant in (Q). The quantities l and 
£ in (^) also need to be zeta- regularised - this was done in p!l| , where it was 
shown that the regularised £ is given by 

C = dimff°(M) - dimH 1 (M). 

So in the present case, where H°(M) = R and H 1 (M) = 0, we have 

C = 1-0=1. (8) 

Using the formulae d\ ~ *di* and ** — 1 (modulo a possible sign), we get 
di = *d\* and therefore (*di) 2 = *di * d± = d*di, which gives 

det'((*di) 2 ) = det'Kdi). (9) 

Substituting (|) and (|) in (|) we get 

Z{\) = e-T'(_)-!det'(cM^det'(dScZo)* det(#5fo) - *. (10) 
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We now rewrite the product of determinants in (|10|) in terms of the Ray-Singer 
torsion p(]] of M. Since the Hodge star operator * is unitary with ** = 1 and 
c?o = *d,2* (modulo a possible sign) we have 

det'(do^o) = det'(*d 2 * do) — det'(*(*d 2 * do)*), 
= det'(d 2 *d *) = det'(d 2 d 2 ), 
= det' (d* 2 d 2 ). 

It follows that 

det'(d 'do)3det'(d* 1 di)-i = (dct'(d*d Q )? det(d5di)-*det'(d5da)*)*. (11) 

It is possible to rewrite det(</>Q0o) using a standard result of manifold theory in 
a different form(see Ell]). We start by noting that the integration map 



H 3 (M) 



is an isomorphism, i.e. for each r £ K there is a unique class [a] £ H 3 (M) such 
that J a — r. (Note that the integration map is well defined on H 3 (M) since 
Im a + dfi — f M a i.e. J M d/3 = by Stokes theorem.) Also from the definition 

H 3 (M) = n 3 (M)/lm(d 2 ) = (Im(d 2 ) ® Im(d 2 ) ± ) /lm(d 2 ) = Imfa^, 
it follows that the map given by 

Im(d 2 ) X — R (12) 
is also an isomorphism. Now define the map 

03 : R -> lm(d 2 ) ± 



to be the inverse of (12). Then using the properties of the Hodge star operator 
it can be shown that 

det(<p* 3 fo) = det^or 1 - 

It follows that 

det(#to)-* = (det^ofo)-* detfoS^)')*. (13) 
Substituting @ and @ in @ we get 

Z(A) |= (-)-5r fl s(M)l, (14) 

7T 

where 

tbs(M) = det(0*0 o )-5 det(<?!)*< ? !)3)5det'(d*d )5det'(^di)-5det'(d*d 2 ). (15) 
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This quantity trs{M) is the Ray-Singer torsion of M [ p0[ . It is a topological 
invariant of M i.e. is independent of the metric of M. Thus the modulus | Z(X) 
of the partition function, given by (p^)-(pl|), is a topological invariant. 

We are now ready to construct a discrete version of the preceding topolog- 
ical field theory which reproduces the continuum expression for the partition 
function where subdivision invariance is the discrete property corresponding to 
topological invariance. We will see that in order to do this it is crucial that 
there is an analogue of the Hodge star operator in the discrete theory. As we 
will see in the next section, this requires a field doubling. Therefore we consider 
a doubled version of the preceding theory, with the fields uj\ and u>2 hi fi 1 (M) 
and with the action functional (||) changed by 

S(w) = A < w, (*di)w >^ SK wa) = A < (£), C dl t)C 2 ) > ■ ( 16 ) 

The reason for this specific choice of action S(lji , 0J2) for the doubled theory will 
become clear in the next section. An obvious generalisation of the proceeding, 
with T — » T = (^q 1 ) shows that the partition function of the doubled theory 



Jn~ L (M)xQ 1 (M) 



ia l (M)xn l (M) 
can be evaluated to obtain the square of ( |l5| ) 

Z(X) =| Z(X) | 2 = (-)-V fl s(M). (17) 

Note that there is no phase factor here. This is because the quantity 1 = d + +gL 
for the action S in (16) vanishes since T = (J^g 1 ) has a symmetric spectrum. 



4 The discrete version of the topological field 
theory 

We proceed to construct a discrete version of Abelian Chern-Simons gauge the- 
ory. The Whitney map enables the Abelian Chern-Simons theory to be discre- 
tised by replacing the gauge field (1-form) A S fi 1 (M) by the discrete analogue, 
a 1-cochain x S C l (K). 

The most immediate way to do this is to construct the action Sk of the 
discrete theory by 

\S K {x) = XS(W K (x)) = A f dW K (x) A W K (x). 

J M 

This can be shown to coincide with the discrete action for the Abelian Chern- 
Simons theory introduced in ||. This prescription fails however, in the sense 
that the resulting partition function Z^-(A) is not a topological invariant i.e. 
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is not independent of K, and does not reproduce the continuum expression 
for the partition function. We demonstrate this by considering the resolvent 
for Sk obtained in an analogous way to the resolvent of the continuum action 
S described in the previous section. Let Tr- : G n (K) — » ^(K) denote the 
self-adjoint operator on C X (K) determined by 

S K {x)= I dW K (x) A W K (x) =< Tkx, x > . 

J M 

Then 

Ker(T K ) C Ker(df). 
Since for x € Ker(<if") we have 

<T K x,x> = I dW K {x) A W K {x), (18) 

W K (d?x) A W K {x). (19) 

M 

Thus the discrete analogue of the resolvent (||) is a resolvent for Sk ■ 

-> K n°(M) -^-^ Ker(df ) C Ker(T K ) = Ker{S K ) 0. 

The resulting partition function is the discrete analogue of the partition function 
Z(X): 

Z K (X) = det'((^)*0*r W((d*)*<) W(--T K )-i 
In |ll| the following formula for Tk was obtained: 

where the sum is over all 1-simplices [v2, V3] such that [vq, vi,i>2, V3] is a 3-simplex 
with orientation compatible with the orientation of M. 

It is possible to show g that det((^f )*<$) = = dimC (K) = the number 
of vertices of K. Then the failure of the discretisation prescription can be 
demonstrated by showing that the quantity 

I Z{\) | 2 = -Ldet'(af ^)det'(T|)-i 

is not independent of K. 

A discrete version of the doubled topological field theory with action (|l6|) has 
been constructed in []ll| in such a way that the expression ( |l7|) for the continuum 
partition function is reproduced. We briefly describe this in the following. 
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The discretisation prescription is 

(wi.wa) e n\M) x fl x (M) -► 0,y) e C 1 ^) x C 1 ^), (20) 
5(c) = A < (-), C dl t)C 2 ) >- (z,y) = A < (*), („ V o"*)® >, (21) 

where .K" is the simplicial complex triangulating M, K is its dual, C q (K), 
C P (K), d K and d are as described in the previous section. The analogue 
of the Hodge star operator is the duality operator * K . This is a map * K : 
C q (K) C z - q (K) ( and : C P (K) -> C Z ~ P (K) which explains the need for 
field doubling and the expression Q) for the discrete action 5k {x,y)- There is 
a natural choice of resolvent for §K(x,y) : analogous to the resolvent ([|) in the 
continuum case. It is 

The partition function is 

Z K {\) = [ VxVye- §K{xy) . (22) 

Jc 1 (K)xC 1 (K) 

Evaluating this by Schwarz's method with the resolvent above leads to: 

Z K (X) = (£)-i+"o*-"i* det((^)*^)-*det'((^)*df )*det'((df )*df 
det((4)*0*)-idet , ((d^)*d^)4det'((df)*df)-i. 

There is no phase factor in ( p2[ ) since £ vanishes just like in ([H]). We have also 
used the fact that £ = 1 — Nq + iV/^ which is shown in [[llj . 

Now rewrite the determinant involving K - objects in terms of determinants 
of A'-objects. Modulo a possible sign ± we have the formulae [flTf 

(**)-* = (**)* = **, (23) 

(* K )- X = (* K )* = **, (24) 

«)* =*^- 9 -i* K , (25) 

(<)* =**-d£_p-i**. (26) 

( The ± signs are omitted because they will all cancel out in the following 
calculation.) Now: 



det'((d^)*df) = dct'(* K df **df), (27) 

= det'(**(* K df d^)*^) (28) 

= det'(df * k * K ), (29) 

= det'(df (df )*) (30) 

= det((df )*df ), (31) 
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and 



det'((df)*df) = det'(* K < **df), (32) 

= det'(**(**df *^ df )* K ) (33) 

= det'« * x ), (34) 

= det'(df (df )*) (35) 

= det((df )*df ). (36) 

The integration map (|l2|) has a discrete analogue 



Ker(df)- 1 -> R, 

a -v <o, [M]>, (37) 

where [M] s C${K), the orientation cycle of M, i.e. the sum of all 3-simplices 
of K, oriented so that their orientations are compatible with the orientation 
of M. ( Note that a S Ker(d^) 1 - C C 3 (K) can be evaluated on any element 
a 6 C${K) to get a real number < a,a >G R.) Define the map 

(j>3 :R—> Ker(df ) ± (38) 

to be the inverse of (|37j). Then using the properties of *^ and * K , it can be 
shown that 

det((0f )*0f ) = det(^)*^)- 1 . (39) 
Now using (|3j]), ( ]3l| ) and ([36]) we can rewrite (p3) as 

^(A) = (-)- 1+ <- N ^r K (M), (40) 

where 

2 

r^(M) =det((^)*^)^det((0f)*^)^ ^ det '(«)*<)^ ( " 1), ■ ( 41 ) 

This quantity tk{M) is the R-torsion of the triangulation K of M. It is a 
combinatorial invariant of M i.e. is independent of the choice of triangulation 

This is the untwisted torsion of M, more generally the torsion can be 
"twisted" by a representation of tti(M). The factors involving the determinants 
det'((d^)*d^) constitute the usual Reidemeister torsion of M When these 
are put together with the factors involving det((0f )0f ), i = 0, 3, as in we 
get the R-torsion "as a function of the cohomology" introduced and shown to 
be triangulation-independent in po|. 
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The expression ( pl| ) for tk(M) is analogous to the expression ( |l5| ) for the 
R-torsion trs(M), and in fact it has been shown pll [2^] that these torsions are 
equal 

t k (M) =t rs (M). 

It follows that the partition function (|4^) of the discrete theory coincides with 
the partition function (|l^) of the continuum theory, except for the Zf-dependent 
quantities Nq and appearing in (|4^) . These quantities can be removed by 
a suitable if-dependent renormalisation of the coupling parameter A. 
It is possible to show that Q 

det((<^)*^) =N« 
det((0f)*0f) = 

We will use this result in our numerical work. 



5 Numerical Results 

We are now in a position to proceed to numerically evaluate the discrete expres- 
sions for the torsion obtained. This allows us to check the underlying theoretical 
ideas by numerically verifying that the discrete expressions agree with expected 
analytic results. It also allows us to check that the results obtained are sub- 
division invariant. The subdivision invariance of torsion is demonstrated by 
showing that if any simplex of the triangulation is subdivided, the value of the 
torsion does not change. This is what is meant by topological invariance in the 
discrete setting. The expected analytic result for torsion for a lens space L(p, 1) 
is (See @) 

T(L(p,l)) = i. 

Thus T(S 3 ) = T(L(1,1)) = 1. 

We also show, numerically, that the discrete expression for the Chern-Simons 
partition function obtained without using the * operator is not a topological 
invariant. This shows very clearly the importance of the doubling construction 
method used in the discretisation method, for capturing topological information. 

In order to proceed, we need to efficiently triangulate the spaces S 3 and 
L(p,q). First we triangulate S 3 . We do this by considering a four dimensional 
simplex [vq, V\,V2,V3, V4) and observing that the boundary of this object is pre- 
cisely the triangulation K of S 3 which we require. Next we turn to spaces L(p,q) 
which we need to triangulate in order to proceed. An efficient triangulation of 
this space has been constructed by Brehm and Swiatkowski |25) . We use this 
procedure for our computations p6fl. 
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We can now summarise our numerical results. The R-torsion for a simplicial 
complex K, with dual cell complex K, for either S 3 or L(p,l) involves evaluating 



T = \ ~aF~\F det(^do) det^)" 1 det(d 3 d 2 ), 

V ^ v 0^ v 3 

where Ni =nos. of i-simplices in K, 81,82, 83 are boundary operators on K and 
do,d±,d2 are coboundary operators on K. Note that, in the discrete setting, 
these operators can be expressed as matrices. We do this by using the vertices, 
edges, faces and tetrahedra of our complex as a basis list. Any p-simplex in the 
complex can then be expressed in terms of this. Since we know what d maps the 
various basis list elements to, we can set the coefficients of its matrix represen- 
tation. When we say det(c?ido), for example, we simply mean the determinant 
of the matrix which results from multiplying the matrices corresponding to the 
operators 8\ and do- 

If our complex consisted of just one triangle [0, 1,2], the basis list would be: 

[0], 

1 [1], 

2 [2], 

3 [0,1], 

4 [0,2], 

5 [1,2], 

6 [0,1,2]. 

We know that d[0] — [1,0] + [2,0]. This can be expressed in terms of our 
basis list as d acting on basis element going to —3 — 4. So the matrix d for 
this complex is: 



/ 




















\ 


















































-1 


1 



















-1 





1 



















-1 


1 














V 











1 


-1 


1 


/ 



Note that the last column is zero since d has nothing to map a 2-simplex to 
and that the first three rows are zero since nothing is mapped to 0-simplices. If 
we act on with this matrix we get —3 — 4, as expected, and d 2 = 0. 

If we do not use K then from section [|, 

f =| Z(l) | 2 = J F (det'af^)(detTl)-i 

We evaluated the quantity T numerically for various triangulations K of S 3 
and found the following results shown in Table for the change of T under 
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Complex 


Xi 


NEW 


Xi*Xl 
NEW 


N 3 


N a 




5-5 


125 


15625 


1 


5 


5 


0.04 


8-6 


1152 


2985984 


0.444 


8 


6 


0.0123 


9-6 


2304 


15116540 


0.351 


9 


6 


0.00975 



Table 2: Results for T. 



p 


No. of zero-modes of A p 


Dim HP 





1 


1 


1 








2 








3 


1 


1 



Table 3: Zero modes for S 3 . 



subdivision where m — n corresponds to a triangulation of S 3 with n vertices 
and m tetrahedra and where Xi — det 

It is clear that T for S* 3 is not subdivision invariant. We next evaluate T for 
S 3 and for L(p,l) and check that it is indeed subdivision invariant and agrees 
with the analytic calculations for L(p,l), with p=2,3,4 and 5. These results are 
shown in Tables f||(]. 

As a check on the numerical method we also count the number of zero modes 
of the Laplacian operator on the different p-chain spaces. These numbers gives 
the dimension of the Homology groups and are shown in Table ||. 

5.1 S 3 

As a further check, a systematic way of carrying out subdivision known as the 
Alexander moves [^7j was used to study the subdivision invariance of the torsion. 
There are four such moves in 3 dimensions. They are best explained by example. 
We have 

1. [0,1,2,3,4]- > [x, 1,2, 3, 4] + [0, x, 2, 3, 4]. 

2. [0,1,2,3,4]- > [x, 1,2, 3, 4] + [0, x, 2, 3, 4] + [0,l,x,3,4]. 

3. [0,1,2,3,4]- > [x, 1,2, 3, 4] + [0, x, 2, 3, 4] + [0,1, a;, 3, 4] + [0,1, 2,x, 4]. 

4. [0,1,2,3,4]- > [js.1,2,3,4] + [0,a;,2,3,4] + [0, 1, X, 3, 4] + [0, 1, 2, x, 4] + 
[0,1, 2,3,a;]. 
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Complex 


Xx 


x 2 


a 3 


A i * A 3 

X, 


A 3 


N 


T 


s3 


625 


15625 


625 


25 


5 


5 


1 


als3 


5184 


2.985984* 10 6 


2.7648 * 10 4 


48 


8 


6 


1 


a2s3 


7776 


1.511654* 10 7 


1.04976* 10 5 


54 


9 


6 


1 


a3s3 


5184 


2.985984 * 10 6 


2.7648 * 10 4 


48 


8 


6 


1 


a4s3 


625 


15625 


625 


25 


5 


5 


1 



Table 4: R-torsion for S 3 . 



These are all natural operations in that the first (nth) move corresponds to 
adding a vertex splitting the 1-simplex (n-simplex) [0, 1] ([0,1, ... ,n]) and con- 
necting it to all the vertices resulting in two (n+1) tetrahedra. 

The torsion T, in terms of its component determinants, and the way they 
change under the Alexander moves is exhibited in the table below, 
where X; t = det9jC?j_! and a2s3 means the complex which resulted after the 
type 2 Alexander moves were performed on the S 3 . It's clear from this that 

T = -J— det(5 1 d )(det d 2 d 1 y 1 det(d 3 d 2 ) 
JV3JV0 

is subdivision invariant and thus a topological invariant of the manifold. 

As a final check we tried several other subdivisions. We took a given triangu- 
lation and barycentric subdivided one or more (n) of its faces to get ls3, 2s3, . . . 
(ns3). The results are in Tables | and |. 



Complex 


X\ 


x 2 


A 3 


Ai*A 3 

X,. 


A 3 


A 


T 


s3 


625 


1.562* 10 4 


625 


25 


5 


5 


1 


ls3 


5184 


2.985984* 10 6 


2.7648 * 10 4 


48 


8 


6 


1 


2s3 


3.9448 * 10 4 


5.2166* 10 8 


1.018325* 10 6 


77 


11 


7 


1 


3s3 


2.79936 * 10 5 


8.418024* 10 10 


3.36798 * 10 7 


112 


14 


8 


1 


4s3 


1.876833* 10 6 


1.266709* 10 13 


1.032626* 10 9 


153 


17 


9 


1 


5s3 


1.270325* 10 7 


2.037398 * 10 15 


3.207653* 10 10 


200 


20 


10 


1 



Table 5: R-torsion for S 3 . 



5.2 Lens spaces 

We conclude with the results for the lens spaces. The results are shown in 
Table §. As can be seen, these agree extremely well with the known analytic 
result T(L(p,l))=l/p. 
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Complex 


Xi 


x 2 


x 3 


x 2 


N 3 


N 


T 


L(2,l) 


1.062937* 10 1U 


3.618662* 10 29 


3.74484 * 10 21 


110 


40 


11 


1/2 


L(3,l) 


9.108* 10 12 


1.143589* 10 43 


1.08817* 10 32 


86.666 


60 


13 


1/3 


L(4,l) 


1.1027* 10 16 


4.0468 * 10 58 


2.89 * 10 44 


78.7487 


84 


15 


1/4 


L(5,l) 


1.790617* 10 19 


1.761432* 10 76 


7.49187* 10 58 


76.16 


112 


17 


1/5 



Table 6: R-torsion for lens spaces. 

6 Conclusions 

The method of discretisation introduced works extremely well. The main point 
of the method is to construct discrete analogues for the set (fi p , d, A, *). Previous 
work in the direction has neglected the Hodge star operator, * g, We 
have thus demonstrated that the Hodge star operator plays a vital role in the 
construction of topological invariant objects from field theory. We were able to 
construct an expression for the partition function which is correct even as far 
as overall normalisation is concerned. Mathematically the equivalence between 
the Ray-Singer torsion and the combinatorial torsion of Reidcmeistcr was proved 
independently in 1976 by Cheeger |^2| and Miiller [p3| . It is nice to see the result 
emerge in a direct manner by a formal process of discretisation. On the way we 
also had to double the original system so that K, the triangulation , and K, its 
dual, are both present. If this doubling and the reason for it are overlooked, then 
the topological information present in the discretisation is lost, as our numerical 
results demonstrated. 

It is clear that the geometry motivated discretisation method introduced is 
very general and that it can be used to analyse a wide variety of physical systems. 
In the approach outlined we have captured topological features. In applications 
capturing geometrical features of a problem are also very important. We are 
currently investigating this aspect of our approach. 

A limitation of the method is that there is no simple generalisation to deal 
with non-abelian theories. The discrete analogies of d, d^ were linear : there is 
no natural discrete analogue of c?^ :— d + A, with A a Lie algebra valued one 
form. 
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